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We study the structure of the no-core shell model wave functions for Li and C by investigating 
the ground state and first excited state electron scattering charge form factors. In both nuclei, large 
particle-hole {ph) amplitudes in the wave functions appear with the opposite sign to that needed to 
reproduce the shape of the (e, e') form factors, the charge radii, and the B(E2) values for the lowest 
two states. The difference in sign appears to arise mainly from the monopole Ahu = 2 matrix 
elements of the kinetic and potential energy (T-l-V) that transform under the harmonic oscillator 
SU(3) symmetries as (A, /x) = (2, 0). These are difficult to determine self-consistently, but they have 
a strong effect on the structure of the low-lying states and on the giant monopole and quadrupole 
resonances. The Lee-Suzuki transformation, used to account for the restricted nature of the space in 
terms of an effective interaction, introduces large higher-order Ahui — n,n >2, ph amplitudes in the 
wave functions. The latter ph excitations aggravate the disagreement between the experimental and 
predicted (e, e') form factors with increasing model spaces, especially at high momentum transfers. 
For sufficiently large model spaces the situation begins to resolve itself for ^Li, but the convergence 
is slow. A prescription to constrain the ph excitations would likely accelerate convergence of the 
calculations. 

PACS numbers: 21.60.-n,21.60.Cs,21.60.De,24.30.Cz,25.30.Bf,25.30.Dh 



I. INTRODUCTION 

The ab initio no-core shell model (NCSM) permits cal- 
culations of wave functions in very large model-space 
sizes for nuclei at the beginning of the p-shell. For ^Li 
calculations up to IGhu; have been achieved [1]]. Among 
the successes of the model is its predicted energy spectra 
of light nuclei [U, [3, [1] . Towards the end of the p-shell 
a 10?ia; basis calculation has been achieved for A=ll [J. 
In mass 12 the model provides a reasonable description 
of the low-momentum component of the vector and axial 
currents involved in the electro-weak transitions to the 
ground state on ^^N when a three-body interaction is in- 
cluded 0]. In all of these calculations, a Lee-Suzuki Q 
transformation of the nucleon-nucleon interaction is used 
to account for the restricted nature of the space in terms 
of an effective interaction. Group theoretical analyses Q 
of the no-core shell model wave functions have shown that 
the predicted eigenstates of ^^C and ^^O have very large 
overlaps with a small sub-space of the full model space, 
with the sub-space being defined by the most deformed 
symplectic basis states. The purpose of this paper is to 
examine the structure of the umlti-huj terms in the wave 
functions for the low-lying states at the beginning and 
end of the p-shell in more detail. For this we compare 
NCSM predictions with measured elastic and inelastic 
(e, e') charge form factors in ^Li and ^^C. 

The shape of the electron scattering form factors pro- 
vides a direct probe of the magnitude and structure of 
the higher shell components in the wave functions. The 
charge form factors have the additional advantage that 



two-body meson-exchange currents do not contribute sig- 
nificantly to the form factors below about 2 fm"-'^ [Sj. In 
all (e, e') calculations presented here we use a bare one- 
body operator; as discussed below, the introduction of 
an effective operator (to compensate for the truncated 
model space) does not significantly affect our conclusions. 

Our first main finding is that the ph amplitudes that 
contribute significantly to electron scattering appear with 
the opposite sign to that needed to replicate the experi- 
mental form factors, elastic and inelastic, and the charge 
radii. The 2hu} contributions to the inelastic form fac- 
tor change sign, in agreement with experiment, for suffi- 
ciently large model spaces for ^Li; however, higher-order 
terms do not within the model spaces we examined. Sec- 
ond we show that the symplectic (A,/i), Ahu = 2 ph 
amplitudes in the wave functions are sensitive functions 
of the oscillator parameter. 

II. THE ELASTIC CO FORM FACTORS 

The ground state CO form factor is the Fourier trans- 
form of the charge density, and contributions from two- 
body charge operators and/or relativistic corrections are 
negligible for momenta up to about 2 fm^^ [8]. In a har- 
monic oscillator (HO) basis, the Ohio p-shell charge form 
factor is given by 

Fop-Op{q^) - \/3(l - 2/iy)exp(-y) (1) 

where y = {bq/2)'^ and b is the oscillator parameter. 

When additional shells are added to the model space 
the new contributions to the form factor fall into two 
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main classes. The first of these are the in-shell con- 
tributions (e.g. IsOd — IsOd, IpOf — IpOf, etc.) deter- 
mined by the occupation numbers for the higher shells, 
and the second are from cross-shell ph excitations (e.g. 
Os — Is, Op — Ip, Os — 2s, etc.). At low q the form factor 
is determined by the charge radius 

F{q') = 1 - ll!^ + 0{q') (2) 

and the higher in-shell contributions can be shown al- 
ways to add constructively to the charge radius. The 
cross-shell excitations (which for HO wave functions con- 
tribute to the charge radius only for Ahuj — 2 ph exci- 
tations across two shells) can add constructively or de- 
structively. In a 2hijj calculation for p-shell nuclei there 
are two possible cross-shell contributions, namely, 

Fos-isiq'') ^ VV^y eM~y) (3) 



Fop-ipiq^) - v/l073y(l - 2/5 y) exp{-~y) 

Both of these transform under SU(3) as (A, (x) — (2, 0) 
and represent the 2hu) symplectic contributions to the 
form factors. When these and/or higher shell ATiw = 
2 ph excitations appear with a sign so as to enhance 
the predicted charge radius, they pull in the charge form 
factor in momentum space. 

A. CO form factor for "Li 

Both elastic and inelastic scattering from ^Li have been 
studied extensively [1, [l^, [ll|, [H, [H, ■ Our calcula- 
tions use the CD Bonn nucleon-nucleon interaction p^ . 
The predicted ground state energy of ^Li is least sensitive 
to the choice of the oscillator parameter h over increasing 
basis size for h = 1.79 fm (^iw = 13 MeV) d, which is 
the value we use in the present (e,e') calculations. We 
note that, because of the increased computational dif- 
ficulties in calculating the transition density matrix el- 
ements needed for the (e,e') form factors, the present 
calculations are restricted to a maximum model space of 
lAhuj, to be compared with the IQhuj model spaces used 
to calculate the energy spectra for ^Li. 

Figure 1 compares the experimental data for the elastic 
form factor for ^Li with the model predictions. For model 
spaces up to \Qhuj, the predicted form factor moves out 
further in momentum space as the basis is increased. For 
the largest model space examined, lAhuj, the tend be- 
gins to reverse. In coordinate space the predicted charge 
density (Figure 2) is enhanced in the interior, with lit- 
tle change to the tail as higher shells are added. These 
trends reflect the structure and sign of ph excitations in- 
troduced as the model space increases. There are two 
issues with the structure of the predicted ground state 
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FIG. 1: (color online)The elastic CO form factor for the 
ground state of ®Li. For all but the largest model space, 
the form factor moves further out in q as the model space is 
increased. Experimental data are taken from jl^ l 
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FIG. 2: (color online) The ground state charge density of "Li. 
As the model space increases the density is enhanced in the 
interior, with minor changes in the tail region. For sufficiently 
large model spaces the strength in the interior starts to be- 
come suppressed with a corresponding build up in the tail 
region. 



ph excitations. First, the A?ia; = 2 ph excitation for 
all shells included add destructively to the ground state 
charge radius. These suggest that the sign of the impor- 
tant symplectic excitations in the wave functions may 
be problematic. In momentum-space these excitations 
pull the form factor out in q. Second, there are large 
Os — > ns excitations which pull the form factor out fur- 
ther at higher q. These effects are shown in Figures 3 
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Point charge radius for ''Li in units of fm 
Model space OHlj 2Tiijj Afiuj IQTiuj lAhcu Expt. 

full model space 2.23 2.08 2.1 2.14 2.21 2.38± 0.1 
ph contributions omitted 2.23 2.26 2.31 2.38 2.41 

TABLE L 



than the constructive interference from the higher in-shell 
excitations. This suggests that the predicted sign of the 
Ahuj = 2 ph excitations that transform under SU(3) as 
(A,/^) = (2,0) inhibits the convergence of the calculations. 



B. CO form factor for ^^C 



FIG. 3: color online) The ground state CO form factor for ''Li. 
The figure displays the effect of particle-hole excitations on 
the predicted form factor by arbitrarily setting the one-body 
density matrix elements to zero. Experimental data are taken 
from [g. 
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FIG. 4: (color online) The ground state CO form factor for 
®Li. The figure displays the effect of Os — > ns ph excitations 
on the predicted form factor. The dashed curves show the ef- 
fect of setting the one-body density matrix elements for these 
excitation to zero (arbitrarily). The Os ns excitations are 
sizable and pull the form factor out further at higher momen- 
tum transfers. 



and 4, where the relevant ph one-body density matrix el- 
ements have been arbitrarily set to zero for the purposes 
of displaying their effect on the shape of the predicted 
form factor. 

In Table 1 we show the contributions to the charge 
radius from in-shell versus cross-shell excitations, and the 
destructive interference from the ph excitations is larger 



The trends seen for the elastic scattering form factor 
for ^Li are also seen in the case of ^^C. Again the calcu- 
lations use the CD Bonn interaction T5]. The Ahu = 2 
ph excitations add destructively to the charge radius and 
pull the elastic CO form factor out in momentum space, 
Figure m We note that the charge radius for ^^C (Table 
HI)) is over-predicted, which in part reflects the chosen os- 
cillator parameter, b = 1.663 fm {huj = 15MeV), which 
minimizes the ground state energy. We will discuss the 
choice and effect of b in section |Vl As in the case of ^Li, 
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FIG. 5; (color online) The elastic CO form factor for the 
ground state of ^^C. As the model space increases the form 
factor moves further out in q, and the dominant higher shell 
terms add destructively to the predicted charge radius. Ex- 
perimental data were taken from [13, |18j, |19|, i2^] . 

ph excitations across more than two shells are large, and 
they act so as to move the form factor out further in q. 
For "'^^C these include both the Os ns and Op np 
excitations. Figure [S] displays the effect of arbitrarily 
setting the ph contributions to the form factor to zero. 



Point charge radius for C in units of fm 



Model space 


Qhuj 2Tujj AfiLu QTiuj Expt. 


full model space 
ph contributions omitted 


3.18 2.99 2.95 2.95 2.32± 0.022 
3.18 3.20 3.25 3.24 



TABLE II: 
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FIG. 6; (color online) The elastic CO form factor for the 
ground state of ^^C. The Ghuj calculations displaying the ef- 
fect of the removal of the ph excitations were obtained by ar- 
bitrarily setting the corresponding one-body density matrix 
elements to zero. 



III. C2 FORM FACTORS 

The C2 form factor is determined by the transition 
charge density; there is no significant contribution from 
two-body meson exchange currents below g w 2 fm~^[l6|. 
We examine the longitudinal form factor for scattering to 
the 2.186 MeV (3+ T = 0) state in ^Li and the 4.44 MeV 
(2+ T = 0) state in ^^C. Data for the ^Li C2 transition 
have been measured to g w 3.5 fm~^ [l3,[H,[i3|- Exten- 
sive data are also available for the 4.44 MeV (2+) state 
in [13, [l, [Sill. 

The most significant contributions to the C2 form fac- 
tors for p-shell nuclei in a (0 -I- 2)?iti; calculation are 
transitions within the p-shell (Op — > Op) and ph exci- 
tations across two shells that correspond to the excita- 
tion of the GQR. The latter transform under SU(3) as 
{X,fj,)L = (2,0)2. For HO wave functions the Op Op 
and GQR form factors are |2l| 

Fp^p{q^) = -^8/Wy exp{-y) (4) 

FcQRiq^) = v/24/15 y (1 - l/3y)exp{^y), 

and as before y — (6g/2)^. If a small admixture of the 
GQR is added to the Qhu state so as to enhance the 
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FIG. 7: (color online) The longitudinal electron scattering 
form factor for the 4.44 MeV (2+ T=0) state in ^^C. As the 
model space increases the form factor moves further out in q, 
in contradiction with experimental and theoretical expecta- 
tions. 

B(E2), the form factor is suppressed at high q. For larger 
model spaces higher powers of y are introduced. 

The shape of the predicted F^ is often displayed in 
terms of the C2 matrix element [2l[. In general, the CA 
matrix element is defined in terms of the form factor F\ 
as, 

and 

C2{q) = B{C2y/^=A + By + Cy^ + ... (6) 

where / = fsNfc.m.exp{-y), fsN is the single-nucleon 
charge form factor [22l and the center of mass correction 
is fern. = exp{y/A) [23]. For a HO basis, the number of 
terms appearing in the polynomial ^ is determined by 
the number of shells included in the calculation. For p- 
shell nuclei the experimental C2{q) matrix elements for 
low- lying states generally decreases with increasing g^, 
i.e., the coefficient ratio A/B < 0, where A> 0. 

A. C2 Form Factor in ^^C 

Our calculations for ^^C include model spaces up to 
6huj. We use the CD Bonn ^ and the AV8' ^ nucleon- 
nucleon interactions, as well as the AV8' plus the Tucson- 
Melbourne TM'(99) 3-body [H] interactions. The oscil- 
lator parameter was taken to be 6 = 1.663 fm (hcu = 15 
MeV). 

Figure [7] shows a comparison between the measured 
and predicted form factors for increasing sizes of the shell 



5 



^ 5- 



U 




Ohco CD Bonn 

2h(n CD Bonn 

-- ehcoCDBonn 
■-■ 4h(nAV8'+TM 



1 



2 

q (fm 



FIG. 8: (color online) The C2 matrix element for the 4.44 
MeV (2+ T=0) state in ^^C. The effect of using a three- 
body 
were 



interaction (AVS'+TM) is smaU. 
taken from [H, [M 1^ 



Experimental data 



model space. At low momentum transfers the calcula- 
tions under-predict the form factor. Above q « 1.5 fm~^ 
the calculations over-predict the form factor, and this 
over-prediction becomes increasingly worse as the size of 
the model space is increased. As the model space is in- 
creased beyond Ohuj, the Op Ip excitations add de- 
structively at low q (and thus destructively to the the 
B(E2) value) and constructively at high q, moving the 
form factor out. The predicted form factor is enhanced 
slightly at small q relative to the Ohu! calculation, and we 
largely attribute this to the Os — + Id excitations, which 
appear with the correct sign. We note this differs from 
the elastic CO form factors where excitations from the Os 
shell appear with the same sign as excitations from the 
Op shell. It should be noted that the model spaces ex- 
amined here are restricted to 6huj and that much larger 
spaces may well show very different trends. 

Figure [5] displays C2(g), which we extracted from the 
measured form factor using b=1.7 fm. The experimen- 
tal C2(g) matrix element steadily drops with increasing 
q^. Our multi-?iijj calculations predict C2{q) to have the 
opposite slope, in large part because of the sign of the 
(Op)~^(lp) excitations in the 0+ and 2+ wave functions. 
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FIG. 9: (color online) The transition charge density for the 
4.44 MeV (2+ T=0) state in ^^C. As the model space is in- 
creased the transition density moves towards smaller r. 



B. The C2 form factor for ''Li 

Figure [TU] displays a comparison of the measured 
and predicted form factors for increasing basis size 
for b = 1.79 fm. These calculations use the CD Bonn 
nucleon-nucleon interaction [15]. At low momentum 
transfers, q < 1.0 fm~^, the calculations under-predict 
the form factor. At g > 1.5 fm"-"^, the larger model spaces 
over-predict the magnitude of the form factor, and this 
over-prediction increases with the model space size. The 
general trends seen with the sign of higher shell contri- 
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Figure [5] displays the corresponding transition charge 
density p(t) for the 0^ 2+ transition. The experimen- 
tally determined p{r) peaks at about 2 fm, while the p{r) 
predicted by the NCSM peaks at about 1.5 fm. As the 
model space is increased, the peak moves towards smaller 
r. 



FIG. 10; (Color online) The charge form factor for the S"*" 
2.186 MeV state using ab initio NCSM CD Bonn wave func- 
tions, where b = 1.79 fm. At this value of b {hut = 13 MeV), 
the ground state energy converges fastest over basis size. Em- 
pirical data are taken from Stanford, Saskatoon, and Mainz 
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FIG. 11: (Color online) The charge form factor is calculated 
with b = 1.94 fm {huj = 11 MeV) for comparison to figure [TOl 
The predictions are suppressed at low q and enhanced at high 
q for all size model spaces. Experimental data are taken from 



B(E2) Values for ^Li in units of e^fm** 


6 (fm) 


Ohiu 2huj 4huj 6huj 8huj lOhu) Uhu Expt. 


1.94 
1.79 


6.84 8.14 8.93 9.89 10.73 11.63 13.4 21.8(4.8) 
4.91 6.25 7.03 8.16 9.14 10.22 12.2 



TABLE III: 



butions in the predicted C2 form factors are very simi- 
lar to those seen for ^^C. However, our ability to go to 
considerably higher shells in the case of ^Li allows us to 
explore these trends in more detail. The peak of the pre- 
dicted form factor occurs at higher q than experiment; 
the biggest shift of the predicted form factor away from 
the observed peak occurs between the Ohco and the 2huj 
model space, and a shift in the peak position to even 
higher q continues until Shuj. Above Sfko the peak po- 
sition of the form factor begins to improve. Above the 
peak, at momenta q > 1.3 fm~^, the additional contri- 
butions from higher shell continue to enhance the form 
factor. For low momentum transfers (below the peak of 
the form factor) we see a slow convergence to a magni- 
tude lower than experiment. This low momentum trend 
is consistent with the trend of predicted B(E2) values, as 
summarized in Table [Till 

We also examined the form factor for a set of calcu- 
lations with a different oscillator parameter, namely, b 
= 1.94 fm {huj = 11 MeV), Figure [IH The b = 1.94 fm 
form factors display similar qualitative behavior as the b 
= 1.79 fm calculations. The peak of the former occurs at 
higher q than experiment and continues to shift outward 
until about lOtuv. For large model spaces the situation 
starts to improve. Above the peak of the form factor the 
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FIG. 12: (Color online) The C2 matrix element with b = 1.79 
fm. To calculate the empirical C2 matrix element (data from 
[i3,[i3,[ii), we used b = 1.70 fm. 

higher shell contributions move the form factor further 
out in q with increasing model space. 

The enhancements for the form factors at large q are 
determined by the sign and magnitude of the higher shell 
contributions in the wave functions. As in the case of ^^C, 
this is most striking in the case of the 2huj configurations, 
where the Op — > IpO/ ph excitations add destructively at 
low q and constructively at high q. The slow convergence 
of the B(E2) to a value smaller than experiment is due 
in large part to the fact that these ph excitations add 
destructively to the matrix element. 

A more detailed understanding of the convergence of 
the lower momentum terms in the form factor with in- 
creasing model space can be seen by examining the C2 
matrix element. To obtain an experimental C2 matrix 
element, we chose an oscillator parameter b = 1.70 fm, 
which is close to the value necessary to give the measured 
rms charge radius. The C2 calculated matrix elements 
are displayed in Figures [T^] and [T31 The experimental 
C2 matrix element decreases with increasing momentum 
transfer in contrast to the predictions of the model. 



We graph the ratio of the coefficients B/A appearing 
in eq. 6 in Figure [Til as a function of basis size. For the 
smaller model spaces this ratio has the wrong sign, but 
as the model space increases the sign eventually changes 
in qualitative agreement with experiment. For b ~ 1.79 
fm, the ratio changes sign between Ahuj and 6huj model 
spaces; for the b = 1.94 fm, the sign switches between 
StiLO and lOhuj. 



Except for the largest model space examined, we see 
an approximate linear relationship between basis size and 
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FIG. 13: (Color online) The C2 matrix element with b = 1.94 
fm. Experimental data are taken from [l^ . [l^ [T^ and a C2 
matrix element extracted nsinp b = 1.70 fm. 
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FIG. 14: (Color online) We fitted the C2 matrix element to 
the polynomial A + By + Cy^ . . .. The ratio of B/A is graphed 
along the j/-axis. Since A is always positive, the change in 
sign is exclusively in B. Experimental data was taken from 
[l2.[l^.[l3| and a B/A ratio was extracted using b = 1.70 fm. 



B/A and C/A, Figs. I14I15I The lihuj model space cal- 
culation suggests that the rate of convergence starts to 



» b - 1 .94 fm 

b - 1 .79 fm 
■ Experiment 




FIG. 15: (Color online) The ratio C/A (eq. (3)) from a fit 
to the C2 matrix element in ^Li as a function of the model 
space. Experimental data are taken from [l^. [l3l. [l4j . and a 
C/A ratio was extracted using b = 1.70 fm. 



FIG. 16: (color online) Comparison of the predicted elastic 
CO form factor for the ground state of ®Li with and without 
the inclusion of a 3-body interaction. 



IV. INCLUSION OF A 3-BODY INTERACTION 

The inclusion of a 3-body interaction leads to an im- 
proved predicted level spectrum in ^Li, particularly for 
the splitting between the ground state and the first 3+ 
state. In addition, the magnetic form factor for the 
0+ —> 1^ transition in -"^^C is significantly improved when 
a 3-body interaction is included This is because the 
predicted form factor is very sensitive to the strength of 
the spin-orbit interaction. However, the present (AL = A 
AS = 0) CA charge form factors are largely insensitive 
to the strength of the spin-orbit interaction and conse- 
quently to the 3-body interaction. In Figure [TBI we com- 
pare the predicted Afuv and 6huj predictions for the elastic 
CO form factor for ^Li. Figure [T7| shows the equivalent 
calculations for the transition C2 form factor to the 3+ 
state. In both cases the inclusion of the 3-body interac- 
tion has little effect on the predicted form factor although 
it does improve the shape at higher q very slightly. Fig- 
ure 7 shows the effect of the 3-body interaction for the 
inelastic C2 form factor of ^^C, which is also very small. 



V. DEPENDENCE ON THE OSCILLATOR 
PARAMETER 

The unexpected sign for some of the higher shell com- 
ponents in our NCSM calculations bears strong resem- 
blance to a similar problem found in standard imilti-huj 
HO shell model calculations. When HO standard shell 
model calculations are extended to include multi-?iw con- 
figurations the lack of self-consistency (in the Hartree- 
Fock sense) causes some of the higher shell c omp onents 
in the wave functions to be unphysical [26l. [27l. [281 [29l. [30l. 
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FIG. 17: (color online) Comparison of the predicted C2 form 
factor for the 3^ state of ''Li with and without the inclusion 
of a 3-body interaction. 
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l3ll . [32, l33| . The main problem arises because matrix ele- 
ments of the kinetic energy (T) and the two-body interac- 
tion (y) across two shells {Ahui = 2) are large and oppo- 
site in sign, and they cannot be calculated reliably. The 
dependence of the C2 matrix element on b implies that 
the magnitude and even the sign of the {ph\T + V\OTilu) 
matrix elements depend on b. These ofF-diagonal matrix 
elements across two shells in turn affect the sign of the 
leading ph excitations in the wave functions, as well as 
all of the similar Ahuj = 2 matrix elements up to the 
maximum shell included in the calculation. 



FIG. 18: (Color online) The charge form factor and the C2 
matrix element for ®Li are calculated for several values of b in 
the 2huj model space. 



calculations. 



B. 



Dependence of the ^^C Form Factor on the 
Oscillator Parameter 



A. Dependence of ''Li Form factors on the 
Oscillator Parameter 

We investigated the effect of the oscillator parameter 
on the predicted form factor within the 2hu} model space, 
using four different values of b ranging from b = 1.66 - 
1.94 fm {hio = 15-11 MeV), as shown in Figure[lHl As b 
varies so does the predicted shape of the form factor: the 
width of the peak becomes narrower, and the peak itself 
shifts to lower momentum values with increasing b. The 
change in slope of the C2 matrix element with b (Figure 
fTSb) suggests that for sufficiently small b, the slope of the 
C2 matrix element will become negative, qualitatively 
agreeing with experiment. But such a small value of b 
would would likely result in very slow convergence of the 



In Figure [19] we display the dependence of the pre- 
dicted C2 matrix element in ^^C on the oscillator param- 
eter. These calculations were restricted to a (0 -I- 2)'hLU 
model space. As the oscillator parameter is varied, 
the value of {ph \ T + V \ QTllo) changes consider- 
ably and eventually changes signs. For sufficiently small 
6(< 1.33/m) the slope of G2{q) becomes negative, in 
qualitative agreement with experiment. 



C. Dependence of the Giant Resonances on the 
Oscillator Parameter 

Hoshino et al. (3^ have pointed out that the problem 
with the sign of matrix elements of (T -|- V) in multi- 
hu) shell model calculations also manifests itself in the 
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q (fm ) 



FIG. 19: (color online) Dependence of the C2 matrix element 
for the 4.44 MeV state in ^'^C on the oscillator parameter for 
the (0 + 2)hui model space. 



predicted excitation of the GQR and GMR. The excita- 
tion energy of the GMR reflects the compressibiHty of 
the nucleus. Both are intrinsic properties of the nucleon- 
nucleon interaction; and, the excitation energy should 
not depend on the properties of the HO well. In Figure 
BUI we show the predicted EO strength for two different 
values of the oscillator parameter for our (0 + 2)Tilu model 
space. The large shift in the predicted excitation of the 
GMR from w 35 McV to w 60 MeV occurs because of the 
change in the off-diagonal (ph \ T + V \ OTiu) matrix ele- 
ments for the two values of the oscillator parameter. An 
analogous problem is seen with the E2 strength, Figure 
[U We note that the sensitivity of the excitation en- 
ergy of the giant resonances to the oscillator parameter 
would likely be considerably less for larger model spaces. 
But our (0 -|- 2)?iu! model space calculations exhibit sim- 
ilar sensitivity to that seen by Hoshino et al. [s^]. As 
in the case of the (e, e') form factors, the problem can 
be directly traced to the Ahuj = 2, (A, ^) = (2, 0) ph 
excitations. 
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FIG. 20: Dependence of the predicted Giant Monopole 
Strength (B(EO)) in ^^C on the oscillator parameter for the 
(0 + 2)7110 model space. 
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FIG. 21: Dependence of the T=0 E2 strength built on the 
ground state of ^^C on the oscillator parameter for the (0 -I- 
2)%^) model space. 



D. Effect of an effective two-body electron 
scattering operator 

Since the model space sizes discussed here are finite, it 
is important to address the issue of the impact of effective 
(e, e') operators and whether they can correct for the sign 
of the ATiw = 2 ph excitations in the wave functions. 
Stetcu et aL^3^ investigated how a two-body effective 
contribution affects the E2 and C2 operators. In the Ohuj 
space, they found that the two-body operator moves the 
form factor in the same direction as the larger (0-|-2)?ia; 
calculations move; that is, the two-body operator results 
in the same unphysical enhancement of the form factor 
at large q. 



VI. CONCLUSION 

We have calculated the elastic CO and the first excited 
state C2 charge form factors and the corresponding C2 
matrix elements in ^Li and ^^C within the NCSM using 
one-body bare operators. These calculations reveal two 
primary findings. First, the magnitude and sign of higher 
shell ph amplitudes in the wavefunction do not behave as 
expected. Higher shell contributions add destructively 
at low q and constructively at high q to the form fac- 
tors, contrary to experimental and theoretical expecta- 
tions. The relative sign of the symplectic (A,/i) = (2,0), 
Ahiv = 2 amplitudes cause them to add destructively to 
the charge radii. The large Os ns (and Op — > np for 
^^C) amplitudes, introduced by the Lee-Suzuki transfor- 
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mation, also affect tlie sliape of the form factors, fur- 
ther increasing the magnitude of the form factors at high 
momentum transfers. In the larger model spaces we ex- 
plored for ^Li, the sign of the 2hL0 contributions to the 
inelastic form factor changes, but convergence onto ex- 
periment is slow. 

The second main finding is the strong dependence of 
the magnitude and sign of the off diagonal /S.Tilj = 2 ma- 
trix elements oiT + V ox\ the oscillator parameter. As a 
result, the observables (the C2 form factor, GMR, GQR) 
also heavily depend on the choice of oscillator parameter. 
Furthermore, the ph configurations in the low lying wave 
functions appear with an unexpected sign. These results 
indicates a lack of self-consistency in the NCSM similar 
to that found in the standard HO shell model. While 
there is no known solution, the effects may be minimized 
by including a Hartree-Fock condition in the calculations 
A prescription [2^, [s^l used in HO shell model calcula- 
tions is to invoke the following condition: 

((n + 2)huj ph\T + V\ nna;)(^'^)=(2,o) ^ q 



Certainly this method bears further investigation. But in 
any case a correct treatment of the symplectic terms in 
the wave functions is crucial to obtaining a realistic de- 
scription of electron scattering form factors within a HO 
shell model basis [35 ]. Invoking eq. (7) and/or another 
prescription to handle the ATiw — 2 and higher ph exci- 
tations may lead the NCSM to have as much success in 
predicting momentum-based observables as in predicting 
energy spectra. 
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